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We study thermal and charge transport in a three-terminal setup consisting of one superconduct- 
ing and two ferromagnetic contacts. We predict that the simultaneous presence of spin filtering and 
of spin- dependent scattering phase shifts at each of the two interfaces will lead to very large nonlo- 
cal thermo- electric effects both in clean and in disordered systems. The symmetries of thermal and 
electric transport coefficients are related to fundamental thermodynamic principles by the Onsager 
reciprocity. Our results show that a nonlocal version of the Onsager relations for thermoelectric 
currents holds in a three-terminal quantum coherent ferromagnet-superconductor heterostructure 
including a spin-dependent crossed Andreev reflection and coherent electron transfer processes. 



Heterostructures of ferromagnets (F) and supercon- 
ductors (S) are presently subject of intense study since 
they show interesting phenomena based on the singlet- 
triplet conversion of pairing amplitudes at the interfaces, 
and the resulting spin-dependent proximity effect. Spec- 
tacular examples are long-range triplet Josephson cur- 
rents due to inhomogeneous magnetic order [1], or due to 
the spin-dependence of the interface reflection and trans- 
mission amplitudes [2], that were confirmed in a set of 
pivotal experiments [3-6] . A multitude of coherence phe- 
nomena are understood in terms of spin-dependent An- 
dreev bound states [2, 7-17], intimately related to spin- 
mixing [18] and spin-filtering effects at interfaces [19]. 

A three-terminal superconductor-ferromagnet proxim- 
ity system also allows to access nonlocal effects. For ex- 
ample, in Fig. 1 incoming electrons (current Ii) can be 
reflected from the interface (/r), or enter the supercon- 
ductor, where each builds a Cooper pair with another 
electron, leaving a hole behind that is retroreflected (so- 
called Andreev reflection). These holes can be transmit- 
ted back through the same interface (Jar), or reflected 
to the other interface, where they are either transmitted 
directly as holes (/car) or as electrons via the same con- 
version process as at the other interface in reversed order 
(^cet) (part of these electrons can also be reflected back 
to the first interface contributing to higher order pro- 
cesses). Nonlocal transport has attracted considerable 
interest due to the latter two processes, called crossed 
Andreev reflection (CAR, electron enters at one terminal 
and hole leaves the other terminal, or vice versa), and 
coherent electron transfer (CET, sometimes called 'elas- 
tic cotunneling', electron enters one terminal and elec- 
tron leaves the other terminal, or the same for holes) 
[20-22]. These processes test the internal structure of 
Cooper pairs, and lead to new interesting physics that 
can be, and has been tested experimentally [23-27]. 

In this Letter we develop a theory for the hitherto 
less explored nonlocal thermal transport in ferromagnet- 
superconductor devices, and show that a nonlocal version 



of Onsager relations [28] holds in both the normal and 
superconducting state. In the superconducting state we 
find a strongly enhanced local thermopower and nonlocal 
Seebeck effect. These effects do not require noncollinear 
inhomogeneities in the ferromagnetic regions or at the 
interfaces (a ubiquitous problem for creating triplet su- 
percurrents [1, 2, 4-6, 29]). Thus, our results should be 
readily observable in experiments and offer a way to ac- 
cess the microscopic spin-dependent parameters. 

In linear response the transport coefficients relating 
charge (energy) currents I q (I £ ) to an applied voltage 
AVj = Vj — Vs or temperature difference ATj = Tj — Ts 
(throughout this Letter j £ {1,2} labels the ferromagnet- 
superconductor contacts, and q = —\e\ is the electronic 
charge) of our three-terminal system are: 
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FIG. 1: (color online) (a) The device consisting of two fer- 
romagnets (to the left and right, blue) and a superconduc- 
tor (in the center, green). Trajectories for electrons (black) 
and holes (red) illustrate possible transport processes in the 
ballistic case, as discussed in the text (white arrows denote 
the spin), (b) Equivalent circuit diagram of the setup shown 
in (a) for the diffusive limit including the coherence leakage 
[41]. The interface parameters are discussed in detail beneath 
Eq. (3). 
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This generalized conduction matrix L contains local 2x2 
blocks in the diagonal, and nonlocal 2x2 blocks in the 
off diagonal. The local and nonlocal thermoelectric co- 
efficients L^J in Eq. (1) give rise to large thermoelectric 
effects in the superconducting state, as we will show be- 
low. In contrast, in the normal state these coefficients are 
typically proportional to the asymmetry of the density 
of states around the chemical potential, which is orders 
of magnitude smaller. Microscopically, spin-dependent 
scattering phases at a ferromagnetic contact produce an 
asymmetry, equal in magnitude and opposite in sign for 
the two spin species, in the superconducting spectrum of 
quasiparticles emerging from the contact. Spin filtering, 
which weights the spin directions differently, can resolve 
these asymmetric components of the spectrum. Both ef- 
fects vanish for spin-independent systems. Consequently, 
this situation is not comparable to the thermoelectric ef- 
fects related to supercurrents discussed in the context 
of normal- metal/superconductor Andreev interferome- 
ters [30-33]. The effects we present persist also in the 
absence of a supercurrent emerging from the supercon- 
ducting terminal. 

We find that the matrix in Eq. (1) (even for non- 
collinear magnetization configurations) is symmetric, 
L = L T , similar to the well-known Onsager symmetries 
[28], however for a nonlocal setup, that contains ferro- 
magnetic leads and includes supercurrents in the super- 
conducting terminal as well as crossed Andreev reflection 
and elastic cotunneling processes between the contacts. 

We begin our theoretical analysis with the descrip- 
tion of the interfaces between the superconductor and 
the ferromagnets. Each conduction channel n between a 
superconductor (S) and a ferromagnet (F) (with homo- 
geneous magnetization throughout the interface region) 
is described by a scattering matrix 

( r e i(p ™ t e iip ™ \ 
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where a G {t, I}, and unitarity requires r^+t^ = 1 and 
(p%v + ( Pna = Pna + Pna modulo 2tt. This leads for exam- 
ple to spin-dependent conductances (spin filtering) char- 
acterized by a polarization V n = (t^ t — t 2 ni ) / {t nj[ + t 2 ni ) 
and a probability for transmission, T n = (t^ + t 2 l ^)/2 
< (1 + \ VnD~ 1 . Concerning the scattering phases, trans- 
port coefficients only depend on the phase shift between 
the reflections of spin-up and spin-down electrons on the 
superconducting sides of the contact, 5(p n = (p% — ip% , 
called spin-mixing angle. Some of the most striking con- 
sequences of the spin-dependent scattering phases are 
triplet pairing [2, 34] or subgap resonances in the noise 
spectral density [11, 35]. Finally, the combination of both 
spin-dependent parameters V n and Scp n leads to thermo- 
electric effects. We use spin-dependent boundary condi- 
tions (SDBC) [2, 7, 12, 36-38] for quasiclassical Green 
functions in the setups shown in Fig. 1. 



Analogously to the spin-independent theory [39-42] 
the system properties in the dirty limit (elastic mean 
free path much shorter than superconducting coherence 
length) are fully described by the isotropic matrix Green 
functions Q c of the contact region [see Fig. 1(6)], and 
Gj (Gs) f° r the ferromagnets (superconductor), that are 
8x8 matrices in Keldysh(g)Nambu(g)spin space. Q c is de- 
termined through a finite element approach, governed by 
a conservation law for matrix currents [41] (see Supple- 
mentary Material for details): J2j %,c +is,c +^Leak = 
with the normalization condition Q 2 = 1. The leakage 
current 2"Leak describes the decoherence of the supercon- 
ducting order parameter due to a finite diffusion time in 
the central region (defining the inverse of the Thouless 
energy £xh)- The spin-dependent matrix currents X^ c 
from contact j into the superconducting contact region 
(denoted c) are obtained from the SDBC. We introduce 
the notation t na = t n + at' n for spin components of the 
transmission quantized along a magnetization direction 
ra. Choosing the spinor basis & = ^j, \E^, — \E^) 
and following the line in Ref. [37] we find to leading or- 
der in t n , t' n , and S(p n a compact form for the SDBC: 

q 2 
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with ij n = tj n + t'- n kj and hj = \®f z ® (rhj&) (f and 
a are Pauli matrices). The tj n and tj n can be related to 

the Tjn and V jn via (t jn +t' jn mj&) 2 = T jn (l + V jn rh j b t ). 
Performing the sums over n, only few parameters remain. 
In terms of the conductance quantum G q = q 2 /h these 
are G 3 = 2G q £ n T jn , Gf R = G q £ n T jn V jn , and G+ = 
2G qJln Sl Pjni as wel1 as ^Th = SThG s /G q . Here, G s is 
the conductance between the contact region and the bulk 
superconductor, fulfilling Ts, c — ^f~[Gs>G c ]- The above 
procedure is correct for 5(fj n ,Tjn <^ 1? covering the full 
range — 1 < Vj n < 1- The equations for Q c are solved 
numerically and the density of states and the currents 
are calculated as function of the parameter set introduced 
above as described in the Supplementary Material. 

In the clean limit (elastic mean free path much longer 
than superconducting coherence length) we apply the 
theory developed in Refs. [10, 12, 43]. In this case, the 
current density at one particular contact can be decom- 
posed into local (depending on the distribution function 
of the ferromagnet at the same contact) and nonlocal 
(depending on the distribution function of the ferromag- 
net at the other contact) contributions: incoming (/j), 
reflected (/r), Andreev reflected (Jar), crossed Andreev 
reflected (/car) and coherent electron transfer /get (see 
Fig. 1). The total current through contact j into the 
superconductor is given by 

tol _ ra ja i ra ra _■_ ra ( a\ 
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with a G {q,s} and contact index j G {1,2}. We con- 
sider two contacts of diameter that are small compared 
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to the superconducting coherence length £q> and to the 
intercontact distance L. Then, quasiclassical trajectories 
connect the two contacts, with contact i seen from con- 
tact j under a solid angle SQj — A\jL 2 ^ where A\ is 
the area of contact i projected on the plane normal to 
the line connecting the two contacts (here, the z axis). 
The current through contact j is proportional to Aj, and 
its nonlocal part is proportional to AlA^/L 2 , as is the 
nonlocal part of the current through contact i. Nonlocal 
contributions enter also /r and Jar, however they are 
the only contributions to /car and /get- Only nonlocal 
contributions, via the trajectory connecting the two con- 
tacts, give rise to thermopower and Seebeck effect in the 
ballistic limit. 

We write nonlocal current contributions as 
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with (S 2 p/Sfl)\ Pl ^ 2 = (S 2 p/Sfl)\ P2 ^ 1 being the differen- 
tial fraction of the Fermi surface of the superconduc- 
tor with Fermi momentum such that the corresponding 
Fermi velocity connects the two contacts, per solid an- 
gle ft. With the deviations of the distribution functions 
from that in the superconductor, for particles 5/ p , and 
holes, 5f h , the contributions to jj = 2 3 \\ ~ 3j,R + jj,AR - 
jj , cet ~h7 j , c ar are e.g. for contact j = 1: ji,i(V) = 2Sf liP , 



Ji,ar(^) 
Ji,cet(^) 
Ji,car(^) 



2|r lt -^^ t r u e^ 7 o7i| 2 ^i,p, (6) 
(*it*u) 2 K| 2 (l7i| 2 + l7o| 2 )«/i,h, (7) 
(^2t) 2 |^i^i2| 2 (l + l7o| 4 r 2 ^ 2 2 ; )^ 2 ,p, (8) 
(ti t ^)>^ 12 | 2 |7o| 2 (r 2 2 t + r 2 n ) 5/ 2>h , (9) 



with 70 (e) = -A/(e + iu), u(e) = VA 2 -e 2 , Tj(e) = 
7or jt r^e^, u 12 (e) = [c - is(e + r 2 A)/cj] _1 , 7l (e) = 
iii 2 [r 2 c + zs(A + r 2 e)/o;], = (1 - TiFi) -1 , with 

c(e) = cosh(u L/ hyp), s(e) = smh.(ojL/hv^). Finally, 
jj(s) in Eq. (5) is obtained by interchanging and 
Scfj — >• — 5^- for both contacts in the expressions above. 
The distribution functions are 



Ak B T s cosh 2 (e/2k B Ts) 



6f jib (-e). (10) 



Equations (4)- (10) are valid for arbitrary transparencies 
and spin polarizations. Nonlocal effects decay when L 
exceeds the scale of the superconducting coherence length 
(£o = hv-p/k-QTc in the clean limit). See Supplementary 
Material for examples. 

The temperature dependence of the superconducting 
pair potential A is taken into account by solving self- 
consistent ly the gap equation in weak coupling BCS the- 
ory (with its zero temperature value denoted Ao). 

As shown in the Supplementary Material, in ballistic 
systems only processes that involve the opposite contact 
contribute to the local thermoelectric coefficients L^J 



and LjJ. The term 7" AR does not contribute because 

ji,AR(— s) cancels the corresponding term for j'i,ar(^) in 
the expressions for the thermoelectric coefficients [both 
have the same pre-factor (£1^14) 2 , i.e. spin filtering is 
not active here]. In contrast, the expression for I? R does 
not show such a cancellation when contact 1 is spin polar- 
ized, due to the asymmetric combination of transmission 
and reflection coefficients in Ji,r(s) (i.e. spin filtering is 
active) and the presence of spin mixing (5(fi). It does, 
however, require in addition that r 2 ^r 2 ^e^ 2 7^ 1 (which 
means the presence of a second contact) in order for it to 
cause nonzero thermoelectric effects. When the impurity 
mean free path or the dimension of the superconduct- 
ing terminal shrinks below direct backscattering due 
to impurities or surfaces contributes and leads to a local 
thermopower even in a two- terminal device. 

As the mechanism behind the thermoelectric effects 
can be understood from the density of states (DOS) in 
the contact region, we discuss first this quantity. In the 
dirty limit (see Fig. 2) for = the DOS displays 
peaks at e = A resulting from the superconducting leads 
and the proximity induced minigap. The magnetization 
directions are chosen parallel. Increasing simultane- 
ously in both terminals leads to a Zeeman splitting of 
the minigap in spin-up and down parts and consequently 
breaks the symmetry of the spin-projected DOS (SDOS) 
around the Fermi energy (see Fig. 2). Hence, we ex- 
pect a nonvanishing thermopower if a spin-filtering term 
qMR - g p resen t simultaneously. An equivalent discussion 
of the SDOS depending on the spin-mixing angle Sep for a 
ballistic system is done in [13]. The subgap peaks there 
are much sharper compared to the washed-out peak in 
the dirty limit. This is associated with the fact that only 
trajectories connecting the two contacts contribute to the 
nonlocal transport, in which case it is governed by a sin- 
gle length L. This is not the case in diffusive structures, 
where quasiparticles take random paths of various length 
between the contacts (and back to the same contact). 
Nevertheless, both ways lead to an asymmetry in the 
SDOS and consequently to the astonishing prediction of 
giant thermoelectric effects for spin-polarized interfaces. 

We now turn to the experimentally relevant question 
how to define a nonlocal thermopower S12 = — AVi/AT 2 , 
which is not unique in contrast to the local thermopower 
Sj = —AVj/ATj = L^J/(T S L]J). In the Supplementary 
Material we discuss several possibilities to relate voltage 
and temperature differences between the two ferromag- 
nets and the superconductor avoiding a control of en- 
ergy currents. In this Letter we chose to define the ther- 
mopower at contact 1 via S12 = /(T$Ly[), which 
is caused by a temperature difference AT 2 at contact 2 
under the conditions AV 2 = 0, ATi = 0, If = 0. 

In Fig. 3 we show the dependence of S = <Si 2 on the po- 
larization and spin mixing for T/T c <C 1 assuming equal 
ferromagnets. The clean and the diffusive limit show 
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FIG. 2: (color online) Density of states D in the contact re- 
gion for d = G 2 = O.IGs, R = Gf R = 0.005G S (10% po- 
larization) and r]Th = £ThGs/G q = 0.5Ao (with the Thouless 
energy €Th of the contact region), (a) Total DOS depending 
on the spin-mixing term G^ for equal ferromagnets. The G^ 
term splits the pseudo gap into the different spin directions. 
(b) shows the asymmetry in the SDOS for spin-down (the 
spin-up SDOS looks equal but mirrored at the e = axis). 



similar behavior, in particular for weak polarizations. 
For large polarization, values of more than 100/iV/K are 
achievable in both limits. Both limits exhibit the same 
point symmetry with respect to the origin, and vanish 
if one of the spin-dependent parameters vanishes. This 
behavior is understood from the SDOS as follows. The 
symmetry of S with respect to the origin is according to 
Eq. (3) a consequence of a tt rotation in spin space. The 
trace in the current formula (shown in the Supplementary 
Material) is invariant under such a unitary transforma- 
tion. The sign change with respect to the axes can be 
understood by Fig. 2. The two spin projections produce 
thermoelectric effects with opposite signs. Depending on 
positive or negative G MR one or the other of the two con- 
tributions will be weighted more. Thus, a sign change 
in G MR changes the sign of the thermopower. On the 
other hand, a sign change in interchanges the roles 
of spin-up and spin-down contributions to the DOS, and 
hence changes the sign of the thermopower too. Similar 
arguments explain the zero crossing of the thermopower 
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FIG. 3: (color online) Nonlocal thermopower S = 
/TsLlY for a symmetric setup as function of polariza- 
tion V and spin-mixing parameter in the clean (a) and the 
dirty (b) limit for T = Ts = 0.1 T c . We assume equally po- 
larized channels, V n = V. In (a) T n \ = 7i=0.1 = 72 = T n2 , 
L = 0.5^o, SCli = Sn 2 = tt/20; in (b) Gi = G 2 = O.IGs 
and r]Th = 0.5A . S is plotted in units of gk^/\q\ where g = 
-T 2 (1+V 2 )SQ 2 /2tt in the clean limit, and g = -G 2 /{G 2 +G S ) 
in the dirty limit. 
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FIG. 4: (color online) Temperature dependence of local and 
nonlocal thermoelectric coefficients for a symmetric setup in 
the clean and dirty limit for various spin mixing parameters 
Sip and G^ . Both coefficients are normalized to the normal 
state value of the nonlocal conductance (L^)t>t c , and are 
plotted in units of kBT c /\q\. Here, V n = V — 0.1, ?7xh = Ao, 
and all other parameters are the same as in Fig. (3) 



when both spin-polarized peak positions in Fig. 2(a) cross 
the Fermi level. The same mechanism leads to a sign 
change in the clean limit, when the spin-split Andreev 
levels cross at the Fermi energy. Here the effect is even 
more drastic since the width of the crossing peaks is de- 
termined solely by the transmission to the ferromagnets. 

We determine the coefficient matrix L in Eq. (1) for 
temperatures across T c . We concentrate on the param- 
eters Lf[ and L^, as they are representative for local 
and nonlocal thermoelectric properties. In Fig. 4 we 
plot these parameters for different spin-mixing angles and 
10% polarization. Remarkably, we obtain qualitatively 
comparable behaviors of both limits although they are 
based on very different assumptions. The quantitative 
differences are related to the different shifting mecha- 
nisms of the subgap peaks already pointed out above. 
Hence, the best comparison is found for small values of 
Sep (ballistic) and G^ (diffusive). We find a zero crossing; 
at a finite temperature in both cases. The similarity of 
local and nonlocal parameters for small temperatures can 
be understand from the thermally insulating behavior of 
superconductors at small temperatures. 

We observe that the coefficients in Eq. (1) fulfill a gen- 
eralized Onsager symmetry. Onsager's symmetry for lo- 
cal currents was originally derived from microscopic re- 
versibility [28]. Generalizations of Onsager's reciprocity 
theorem have been recently discussed using statistical ar- 
guments [44-46] . Here we find a generalization for nonlo- 
cal superconductor/ferromagnet three-terminal devices, 
that include supercurrents as well as crossed Andreev 
reflection processes. This follows directly from the ana- 
lytical formulas (5)- (10) in the clean limit using relations 
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like Tj(s\— 5(p) = — e;5ip) (an example is given in 
the Supplementary Material), and is verified numerically 
also for the diffusive case. This Onsager symmetry holds 
for any relative angle between the magnetization axes of 
the two ferromagnets. 

In conclusion, we have opened a way of utilizing ther- 
moelectric effects in superconducting spintronics. This 
possibility of controlling energy flow in superconducting 
heterostructures with spin polarized electrodes allows of 
a multitude of novel applications. Particularly interest- 
ing for applications is our finding of a zero crossing in 
the Seebeck coefficients as function of temperature, spin 
polarization, and the relative angle of the magnetization 
axes. This not only would give a possibility to measure 
spin-filtering parameters and the experimentally so far 
inaccessible spin- mixing parameters, but would also al- 
low for sensitive and controllable thermal elements in su- 
perconducting circuits. 
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SPIN DEPENDENT QUASICLASSICAL THEORY 
IN THE DIRTY LIMIT 



In the following we will give a more detailed descrip- 
tion of our calculations. In the stationary case the non- 
equilibrium Keldysh Green function in Fourier presen- 
tation g(r, f, e) = / dt/hG(r,r',t-t')exp (ie(t - t')/H) 
reads: 



G R (r,r',e) g K (r,r',s) 
G A (f,f',e) 



with 



$ R/A (r,r',t-t') = TiO(Ht-t'))({*(t^,y\t\r')}} 
g K (r, r',t-t') = -i( f), &(t\ r ')] ) . 

Here R/A/K labels the retarded, advanced, and Keldysh 
part respectively. We choose the spinor basis 



t(*, r) = (*{ (t, f) s *t (tj ^ ^ f)> _ $t(tj ^) . 



In the following, a and Jo denote the vector of Pauli 
matrices and the unit matrix in spin space, respectively, 
and f z denotes the third Pauli matrix in Nambu space. 
In this basis the quasiclassical isotropic Green functions 
of a bulk BCS-superconductor (S) and a ferromagnet (F) 
are given by (5 > 0, S — » 0) 



g s 



±sign(e) 



(e ± iS) 



^(e±i5)-\A\ 2 V - A 



) cr 



gR/A = ^ ^ 

Here, A is the superconducting order parameter. Note 
that the definition of Of is equal to the one of a normal 
metal. Ferromagnetism enters only in the SDBC via spin- 
dependent phase shifts and interface polarization effects. 
According to, e.g., review [1] the Keldysh component is 

g K = g R h-hg A , 

with the distribution matrix 

Aanh 

' 2T 



h = 



2T 

tanh £ -±^ 



) <T . 



Following the steps done in Ref. [2] the SDBC in the tun- 
neling limit follow only for the limit of small V n . Thus, 
we derived a version of the boundary conditions following 



[3], which is valid for general V n . In linear order in T n 
and Scp n we obtain Eq. (3) of the Letter, which is written 
explicitely as [k = f z <S> {frier), where the ferromagnet is 
described by its magnetization direction unit vector rh] 

2i L (e) = 2G q J2 n [i L n R G R {e){i L n R y - i5<p L n K,Q L {e)] 

=G° [g R (e),g L (e)] +G MR [{k, G R (e)} , G L (e)] 

+ G 1 [f£ R {e)k,$ L {e)]-iG* [k,G L (s)] 

q 2 /h is the conductance quantum) with 

G° 



(G 



G 1 

qMR 

G 4> 



G 9 E„T„(l 
G q E n Tn (l 



Or/l are the Green functions of the right (R) and left (L) 
side respectively. Note that as g^ A ^ K commutes with ft, 
the terms with G° and G 1 can in our case be combined 
to one term of the form G [Of, c ] with G = G° + G 1 = 
2G q ^2 n Tn, leading to a simplified boundary condition 
for F/S contact j [for notational simplicity we use for 
each F terminal the same index as for the corresponding 
contact, j, e.g. we write Gj(e) = QpCVj. 



Gf*{R 3 ,g 3 is)}-iG*k 3 Me) 



For the contact to the bulk superconductor the bound- 
ary condition is is,c(e) = \Gs [Gs(z), Gc(e)] • Dephasing 
in the contact region is described by a leakage terminal, 
i L eak(s) = (G q /4£ T h)[GL^(e),g c (e)], with g^(e) = 
—i(e ± iS)f z (Jo and Thouless energy £xh- The Green 
function of the contact region, Q c , fulfills a Kirchhoff rule 
[4] 2Leak + is,c + ]Cj ^j,c — 0? an d hence is is determined 
from equations of the form [.M,£f c ] = with normaliza- 
tion condition g 2 = 1. The charge and energy currents 
are obtained from the Keldysh component of the matrix 
current [5]: 

I1 = ^jTv[{f z ®a )if c {e)]de 

The density of states in the contact region is defined from 
its retarded Green function like 

D(e)= l -Tr[f z ®<j g*{e)}. 
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CLASSIFICATION OF LOCAL AND NONLOCAL PROCESSES FOR BALLISTIC THREE-TERMINAL 

SYSTEMS 

In this section we show examples for transport processes that contribute to the various terms of Eqs. (6)- (9) of 
the Letter. These terms involve propagation of particles or holes, which will be represented as full lines and dashed 
lines in the figures following. Certain processes involve a conversion between particles and holes, accompanied by the 
creation or destruction of a Cooper pair. These processes will be represented by a loop, where a full line turns into a 
dashed line or vice versa. These loops in the following figures correspond to the factors |7o| 2 in Eqs. (6)-(9) of the 
Letter. Propagation between the left and right interface is represented in these equations by the factors U12 and i^i- 
Vertex corrections v\ and V2 correspond to multiple Andreev reflections possible at both interfaces. The factors 71 
and 72 combine propagation between the two interfaces with Andreev reflections at the other interface. 




FIG. 1: Various contributions to the reflection components ii,R. These processes are characterized by effectively scattering a 
particle into a particle or a hole into a hole at the same interface. Only the leading terms are shown, with up to two Andreev 
reflections [denoted as the loops turning particles (full lines) into holes (dashed lines) or vice versa]. 




FIG. 2: Various contributions to the Andreev reflection components 7i,ar- These processes are characterized by effectively 
scattering a particle into a hole or a hole into a particle at the same interface. Only the leading terms are shown, with one 
Andreev reflection. 




FIG. 3: Various contributions to the coherent electron transfer components 7i,cet- These processes are characterized by 
effectively scattering a particle from one interface into a particle at the other interface or a hole at one interface into a hole at 
the other interface. Only the leading terms are shown, with up to two Andreev reflections. 




FIG. 4: Various contributions to the crossed Andreev reflection components /i,car- These processes are characterized by 
effectively scattering a particle from one interface into a hole at the other interface or a hole at one interface into a particle at 
the other interface. Only the leading terms are shown, with Andreev one reflection. 



We underline that Eqs. (6)- (10) of the Letter were obtained by calculating the full expression for the current through 
each interface. This expression can then be split into the various terms in a natural way, with the interpretation as 
shown in the figures above. 
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ANALYTIC PROOF OF ONSAGER'S 
SYMMETRY IN THE CLEAN LIMIT 

Concerning the Onsager symmetry pointed out in the 
Letter we will show in the following how the symmetry 
follows analytically from the formulas (5)- (10) of the Let- 
ter in the clean case. As an example we will choose the 
nonlocal relation — L 2 Y . Using 



dfi, 



dfi,p 



s Q = dfi,h 

ij Ak B T s cosh 2 (e/2k B T s ) d[AV 



dfi,h 



Ak B T s cosh 2 (e/2k B T s ) ~ d[^] ' 



it is seen that in this case the contributing terms to the 
current are the coherent electron transfer and the crossed 
Andreev reflection part, i.e. the part — I^qet + ^car °f 
the total current in Eq. (4) of the Letter. In the following 
we factor out the prefactor f3 = ^nytL 2 fr° m Eq. (5) 
of the Letter. In a first step we write 



T qT _ jqT 
-1^12 — L io 



teV 
^21 



12, CET 



' ^12, CAR 



/3 



/ 



- [+^4l2,CET(»] + [-^12,CAR(^)] , 

a 7) as 

Ak B T s cosh 2 (e/2k B T s ) 



teV , teV 

^21, CET ' ^21, CAR 



I 



~ [+g^21,CET(g)] + [+^21,CAR(g)] ^ 

4k B T s cosh 2 (e/2k B T s ) 



where we introduced A jiiCET/CAR = (jj,cET/CAR + 
jj,CET/CAR)/^/i, P /h following the definition of the cur- 
rent. Let us first compare the terms for the CET contri- 
butions. The Onsager symmetry for the CET contribu- 
tion follows from Ai2,cet(^) = ^-21,cet(^)- Therefore it 
suffices to prove this latter equality. The only terms in 
the expressions for Ai2,cet and A2i,cet (see Eq. (8) of 
the Letter) that change under interchanging 1 o 2 are 
the terms |^i^i2 1 ? which change into \v 2 u 2 i 1- Following 
the definitions of the Letter, a simple algebra results into 



We find now from the definitions of v\ and u\ 2 and 
the just mentioned relations (we use the abbreviations 

Stpi = {5(p u S(p 2 }, Sipi = {-S(p u -S(p 2 }) 

\viu 12 \ 2 (e;-S(pi) = \vi(e; -5(pi)u 12 (e; -5(pi)\ 2 
= \vl(-e;5(pi)ul 2 (-e;5(pi)\ 2 

= |vi(-e;<fyi)ui 2 (-£;<fyi)| 2 
= |^i^i 2 | 2 (-^;^i) 

From the considerations above for the CET contribution, 
it follows that |^i^i2 1 2 (^; Sipi) = |^2^2i| 2 (^; Sipi). Fur- 
thermore, |7o(e)| 2 = \jo(—e)\ 2 . In order to show the 
relation Ai 2 ,car(^) = ^21,car(-^) we use Eq. (9) of 
the Letter. By introducing = (ti^i) 2 (r 2 ,^ + r 2 ^) 
and B it = (t u t 2 ^) 2 (r| ; + r 2 t ), we find 

^4i2,car(^) = 

= 5 n |vi^i27o| 2 (^; S<pi) + £ it K?ii27o| 2 (£; -fyi) 
= 5 n |^i^i27o| 2 (^; Sipi) + £|tl^i27o| 2 (-£; Sipi) 

and 

^21,CAR(^) = 

= B it \v 2 u 21 <y \ 2 (e; 5<pi) + B n \v 2 u 21 >y \ 2 (e; -5<pi) 

= ^ t |vi^i27o| 2 (^; Sipi) + ^nl v i^i27o| 2 (-^; Sipi) 

and comparing the last lines of these two sets of equa- 
tions, it follows that 

^21 5 CAR(— £) = ^12,CAr(^)- 



V1U12 



= {c(l - 



rir 2 ) - i-[e(i + rir 2 ) + A(r x + r 2 



This expression is the same as the one resulting from 
^2^21? as is also obvious from its symmetry with respect 
to interchanges 1 o 2. Consequently, |^i ^12 1 = 1^2^21 1, 
hence Ai2,cet = ^21,cet? and therefore 



L 



qT 

12, CET 



L 



eV 

21, CET- 



Slightly more involved is the proof for the CAR contri- 
butions. We will prove that Ai 2 ,car(^) = ^21,car(— £)• 
We will use the relations Tj(e; —5(pj) = —Tj(—e;Sifj) 
(the semicolon separates arguments from parameters), 
which follows from 70(e) = — 7q(— e) and the definition 
of the Tj(s;S(pj). Note further that from its definition 
it follows that Tj obviously doesn't change under t^i- 



In a last step we find 



T qT 

^12, CAR 



~gg^l2,CAR(g) 

4k B Tscosh 2 (e/2k B Ts) 

+°° -ggA 21 , C AR(-g) 

00 4fc B T s cosh 2 (e/2k B T s ) 
+ °° e'qA 21t cAR(e') 



ds 



ds 



OO 



4k B T s cosh z (e'/2k B T s ) 



ae — 1^21, car • 



All remaining Onsager relations are proved analogously 
Y)j(s\ —Sipi) = — 7j(— S]S(fi) will be needed, and for the 
Andreev contributions Jar the symmetry requires inte- 
gration over £, similar as for /car above]. 

Note that for symmetric contact parameters B^ = 
B n , and consequently Ai 2 ,car(^) and A 2 i 1 car(s) are 
both even functions in e, leading to L^ CAR = ^I^car = 
0. For antiparallel alignment of the magnetization of 
the two contacts, and otherwise identical contact param- 
eters, T 2 (e;5(p 2 ) = Ti(e\—Scpi) = — r*(— s; <fy?i), and 
\viUi 2 \(e) = |vi^i2 1 (s), which together with |7o(^)| = 



reV 



^21 — ^21 



0. 



|7o(— s)| shows that L\ 2 = L\ 

From the above symmetry relations it follows that for 
identical contacts and parallel magnetization, there are 
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only 3 local and 4 nonlocal independent coefficients; for 
antiparallel magnetization all nonlocal thermoelectric co- 
efficients are zero, and there are 3 local and 2 nonlocal 
independent coefficients. In the general case, for asym- 
metric contacts, there are 6 independent local and 4 in- 
dependent nonlocal coefficients. 



VANISHING OF PURE ANDREEV REFLECTION 
CONTRIBUTIONS TO THE THERMOELECTRIC 
EFFECTS IN THE CLEAN LIMIT 

Similarly as for the CET and CAR contributions, dis- 
cussed in the last section, for the AR contribution one 
can show that A^ar^) = ^jj,ar(— e). This follows 
directly from Eq. (7) of the Letter, and the fact that 
7o(e) = -7o(- £ ) and 7jM - ^j) = -7j(-e; <%)| 2 : 
introducing B'-^ = (tj^-tji) 2 = Bj we find 

A jjlAR (e) = B 'i,n (I«i| 2 (l7i| 2 + l7o| 2 ) (e;S^) 
+ B' ut (\v 1 \ 2 (\ ll \ 2 + \ 7o \ 2 )(e;-5p i ) 
= ^,n{(K| 2 (l7i| 2 + l7o| 2 )(e;^i) 

+ (M 2 (l7i| 2 + l7o| 2 ) (-£;<%)}, 

leading to Ajj^Aii(— e) — A^ar^)- Introducing this into 
the integrals for Lj? AR and L^ AR shows that both van- 
ish identically: 



rqT 



rsV 



0. 



NONLOCAL ORIGIN OF ALL 
THERMOELECTRIC COEFFICIENTS IN THE 
CLEAN LIMIT 

As shown in the previous section, the Andreev reflec- 
tion contributions to the local thermoelectric effects van- 
ish, such that local thermoelectric coefficients are solely 
due to the /r contribution, determined by Eq. (6) of the 
Letter. The presence of the second spin-polarized con- 
tact is, however, crucial. In the absence of a second spin- 
polarized contact the local thermoelectric coefficients 
etc. vanish, as we show in this section. 

This can be understood in the following way. In the ab- 
sence of a second contact quasiparticles entering through 
contact 1 can either be Andreev reflected, or scattered 
from some surface point of the superconductor. After 
scattering, they will undergo Andreev reflection, with a 
hole being backscattered (we assume that any surface 
point from which direct retro-reflection into the origi- 
nal contact occurs is far away from the first contact on 
a coherence length scale, and thus contributes negligi- 
bly). Typical backscattering events are presented in Fig. 
1 of this Supplementary Material. We show now, that 
backscattering from any specular surface point of the su- 
perconductor will not lead to any thermopower either. 



Assume that the second contact is specularly reflecting 
and has zero spin-mixing angle. We show that the re- 
sulting thermoelectric coefficients vanish. For S(f2 — 
and = = it follows that r 2 (e) = 70(e), which 
leads after some straightforward algebra to 71(e) = 70(e) 
(meaning that the propagator coming from interface 2 is 
now equal to one coming from the bulk in the supercon- 
ductor). Thus, with 



Aii iR (e) = 2 
+ 2 



r lt -r n e iS *i 70(e) 



1 -r n r u e^i 70(e) 2 
r u -r lt e-^ 7o (e) 2 



l-ri t r u e-^7o(e) 2 



(using ri t -^itf t r u e^i7o7i = M r it _r H e '7071]) we 
consider two cases. For \e\ > | A(T) | we have a purely real 
70(e). In this case, 70(e) 2 = 70 (— e) 2 and the expression 
is identical to the one for Ah j r(— e). On the other hand, 
for \e\ < I A(T) | we have |7o| = 1, and can thus write 
70(e) = e**( £ ). As the absolute value does not change 
when multiplying with a phase factor or when taking the 
complex conjugate, in this case we obtain 



Mifrie) = 2 



fit ~ r n e 



70(e) 2 



1 -ri t r u e^7o(e) 2 



1 - ri t r u e^i 7o *(£) 2 



With 7o(e) = — 70 (— e), this expression again is equal to 
^-ii,r( — e). It follows that for all e the relation An^(s) = 
^11, r( — e) holds, and introducing this into the integrals 
for L\i R and Lf^R leads to vanishing results: L\^ R = 
L\Yr = for 5if2 = and t2t = ^24. = 0. 

Thus, nonzero local thermoelectric coefficients (Lf^ , 
L\Y) require the same nonlocal processes as the nonlocal 
coefficients (L^ , LlX)- 



VARIOUS DEFINITIONS OF NONLOCAL 
THERMOPOWER 

In Tab. I we list several possibilities to relate voltage 
and temperature differences between the two ferromag- 
nets and the superconductor avoiding a control of energy 
currents. We assume a temperature difference AT2 is 
applied at contact 2 and a voltage AVi is established at 
contact 1, which determines the thermopower. The su- 
percurrent flowing out of the superconducting terminal 
is given by J| = If + if- I n order to have a clean ef- 
fect, we require no temperature difference at contact 1. 
Furthermore, we avoid scenarios where both voltage and 
current at the same contact should be tuned to zero. We 
also require that either the current through terminal 1 
or through the superconducting terminal should be zero. 
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n 


1 2 


ATi 


AV 2 


ST S 


I 





free 








T QT I j qV 
^12 /^ll 


II 


I 


f + i| = 








r9 T i r <2 T 


III 











free 


j qT T qV T qT T qV 
^12 ^22 _ij 22 ^12 
jqVjqV jqVjqV 
^11 ^22 ^21 ^12 



TABLE I: List of various possibilities to define the nonlocal 
thermopower S = <Si2 = — AV1/AT2 relating the voltage in 
terminal 1 to the temperature difference in terminal 2 in a 
three- terminal device. Control of energy currents is avoided. 

This gives the combinations shown in Tab. I. The exper- 
imentally chosen combination depends on the context. 

In particular, for combinations I-II at each ferromag- 
netic terminal either the voltage or the temperature dif- 
ference is non-zero (but not both), which presents the 
nonlocal effect in a clean way. For case I a finite super- 
current is flowing out of the superconducting terminal, 
for cases II and III there is no such supercurrent. Case 
III is special in that there are no charge currents flowing 
in or out of the device. In case I and II no local ther- 
mopower is present (AV2 = 0). In case III the nonlocal 
thermopower S is accompagnied by a local thermopower 
S 2 = -AV 2 /AT 2 of value 



S^T S 



rqTj-qV 



rqV jqV 
^11 ^22 



jqVjqV 
^21 ^12 



This latter case is interesting from the point of view that 
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FIG. 6: Dependence of thermopower on distance between 
contacts in the ballistic limit. The local and nonlocal ther- 
mopower are plotted for case III. Parameters correspond to 
those in the top row of Fig. 5, i.e. V\ — V2 — 0.99, T1—T2 — 
0.1 and Sif! = 5ip 2 = 0.5tt, <5Hi = SQ 2 = ^(£o/L) 2 . The 
voltage in contact 2 is below that in the superconductor, and 
the one in contact 1 above that in the superconductor. A ther- 
mopower between the ferromagnetic leads is established, with 
no charge currents flowing in the device. The (almost) equal 
magnitudes are a result of symmetric contact parameters; that 
is, differing transmissions lead to differing magnitudes. 



a temperature difference at contact 2 causes thermopow- 
ers at both contact 1 and contact 2 under condition of 
no current flow in the device. These thermopowers ful- 
fill the equation L^S^ in) = L$[(S$ - S^ 2 n) ). In fact, 
numerical calulation shows (see Fig. 6) that the sign of 
Sy^ 1 ^ and S^ 11 ^ is opposite, such that the thermopow- 
ers sum up over the two contacts, thus providing effec- 
tively a ferromagnet/superconductor/ferromagnet junc- 
tion with largely enhanced thermopower mediated by 
spin-polarized Cooper pairs in the superconductor. 

In Fig. 5 we compare results for the thermopower using 
the three cases in Table I. It can be seen that large nonlo- 
cal thermopower arises independent on the magnitude of 
the supercurrent flowing out of the superconducting ter- 
minal. Fig. 6 shows the dependence of the thermopower 
on the contact distance. 



0.2 0.4 0.6 0.8 T/T 0.2 0.4 0.6 0.8 T/T 0.2 0.4 0.6 0.8 T/T 



FIG. 5: Comparison of the three definitions in Table I. In 
the left column in the middle column S^ I:[ \ and in the 

right column S^ 111 ^ are shown. Top row: the spin polarization 
on both contacts are equal and 99% (corresponding to a half 
metal). Bottom row: the spin polarization at contact 2 is only 
50 %: Vi = 0.99 and V 2 = 0.5. For all results Ti = T 2 = 0.1 
and S(fi — S(f2 = 8cp, 8Qi — SQ2 — 7r/20, L — 0.5£o- As can 
be seen, the symmetric case is special for S 1 " 11 ^. Otherwise all 
definitions give qualitatively similar results. 
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